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2004 J. Stat. Mech. P08006] dealing with classical Casimir phenomena in semi-infinite 
wall geometries. In that paper, using microscopic Coulomb systems, the long-ranged 
Casimir force due to thermal fluctuations in conducting walls was shown to be screened 
by the presence of an electrolyte between the walls into some residual short-ranged 
force. Here, we aim to extend the study of the screening (cancellation) phenomena to 
universal Casimir terms appearing in the large-size expansions of the grand potentials 
for microscopic Coulomb systems confined in fully-finite 2D geometries, in particular 
the disc geometry. Two cases are solved exactly: the high-temperature (Debye-Hiickel) 
limit and the Thirring free-fermion point. Similarities and fundamental differences 
between fully-finite and semi-finite geometries are pointed out. 
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1. Introduction 

At zero temperature, fluctuations of the quantum electromagnetic field in vacuum 
manifest themselves via an attraction of two parallel ideal-conductor plates. This 
Casimir effect (for an introduction see pQ) has a universal character in the sense that 
it does not depend on the material constitution of the metallic plates. Casimir's result 
was extended to arbitrary temperatures and general dielectric plates [21 E| > and to ideal- 
conductor walls of arbitrary smooth shapes j3j. For a recent book and review see 0IH]. 
The studied models can be divided according to the geometry of fluctuating walls into 
two basic sets: the semi-infinite systems, in which at least one of the spatial coordinates 
is unconstrained by the walls (e.g. two parallel plates), and the fully-finite systems (e.g. 
a sphere). The applied methods and observed Casimir phenomena usually depend on 
this classification. 

As concerns semi-infinite geometries, in the high-temperature limit defined by the 
validity of the equipartitioning energy law, the Casimir force becomes purely entropic 
[7j; this force is usually called classical since it does not depend on Planck's constant. 
In the purely electrostatic models which do not incorporate the magnetic part of the 
Lorentz force due to charge currents, like the system of scalar photons |H] , the Casimir 
force is divided by a factor 2. 

As concerns fully-finite three-dimensional (3D) conductor systems |H EJ, in the 
high-temperature limit, the Casimir free energy does depend on Planck's constant. 
Furthermore, the presence of both electric and magnetic degrees of freedom is necessary 
for obtaining a Casimir effect. This is no longer true for classical 2D Coulomb fluids 
defined in fully-finite domains. There, the consideration of the pure Coulomb potential, 
defined as the solution of the 2D Poisson equation, leads to a universal Casimir term 
analogous to the one appearing in finite-size expansions of thermodynamic quantities 
for 2D critical systems with short-range interactions among constituents. To be more 
precise, it is known that, according to the principle of conformal invariance, for a finite 
system of characteristic size R, at a critical point, the (dimensionless) free energy has a 
large- R expansion of the form [TU1 HH E21 EE] 



where f3 denotes the inverse temperature. The coefficients A and B of the bulk and 
surface parts are non-universal. The coefficient of the logarithmic Casimir term is 
universal, dependent only on the conformal anomaly number c of the critical theory 
and on the Euler number x of the manifold on which the system is confined. In general, 
X = 2 — 2h — b, where h is the number of handles and b the number of boundaries of the 
manifold (x = 2 for the surface of a sphere, x = 1 f° r a disc, x = f° r an annulus or a 
torus). A simple example is the Gaussian model jI3] which is critical at all temperatures, 
with the conformal anomaly number c = 1. At any temperature of the conducting 
regime, the grand potential of a classical 2D Coulomb system of characteristic size R is 
supposed to exhibit a large- R expansion of type (jl.lj) . with however a + sign in front 



0F = AR 2 + BR - — lni? + const + ■ • • , 
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of the logarithmic term: 

/^Coulomb = AR 2 + BR+^\nR + const + • • • . (1.2) 

o 

Plausible arguments for a critical-like behaviour were first given for Coulomb gases with 
periodic boundary conditions then for Coulomb systems confined to a domain 
by (vacuum) plain hard walls ^H], by inert ideal-conductor walls (i.e. when the 
electric potential obeys Dirichlet boundary conditions) and finally by ideal-dielectric 
boundaries (i.e. when the electric potential obeys Neumann boundary conditions) |18j . 
The explicit checks were done in the Debye-Htickel limit ^H], at the free-fermion point 
of the Thirring representation of the symmetric two-component plasma, based on the 
formalism developed in [20], and for the one-component plasma at (3 = 2 in a disc 
|16[ ITTj . Only recently, a direct derivation of the universal finite-size correction term 
was done, in the whole stability range of temperatures, for the specific cases of the 
symmetric two-component plasma living on the surface of a sphere [21] and in a disc 
surrounded by vacuum [22] . In both cases, the universal prefactor to the lni? correction 
term in (jl.2j) was related to the bulk second moment of the density structure factor 
which is known j2Hj. This is very different from the semi-infinite geometries where the 
universality of the Casimir force is related to the second moment sum rule for the charge 
structure factor [TT1 12~1]. 

It is generally believed that, for semi-infinite systems, the presence of an electrolyte 
between fluctuating conductor walls screens the long-ranged Casimir force to some 
residual short-ranged force [213 f26j. The study of Casimir phenomena via fully 
microscopic Coulomb models has the advantage of a coherent description of electrostatic 
fluctuations inside conducting walls and the image forces acting on the electrolyte 
particles, without any ad hoc separation Ansatz used in usual macroscopic treatments 
[2*5] . We used this strategy to show the screening effect of Casimir forces for semi-infinite 
geometries in paper [21], in what follows referred to as I. 

The aim of the present work is to extend the study of the screening (cancellation) 
phenomena to universal Casimir terms appearing in the large-size expansions of the 
grand potential (jl.2j) for the Coulomb systems confined in fully-finite 2D geometries, 
in particular the disc geometry. The confining disc walls as well as the electrolyte 
inside the disc are modelled by two different microscopic two-component plasmas of 
point-like particles with +/— unit charges in thermal equilibrium. Two cases are solved 
exactly: the high-temperature (Debye-Hiickel) limit (3 — > and the Thirring free-fermion 
point (3 = 2 corresponding to the collapse of positive and negative pairs of point-like 
charges. From the technical point of view, the circular symmetry of the problem leads 
to infinite summations over specific products of modified Bessel functions which have 
to be evaluated by using the asymptotic Debye expansion; we apply a few technicalities 
which help us to simplify a relatively complicated algebra. Similarities and differences 
with respect to screening phenomena in semi-infinite geometries, described in paper I, 
are pointed out. 

The model is defined as follows. We shall consider Coulomb-gas systems of point- 



Screening of universal Casimir terms 



4 




Figure 1. Disc geometry. 



like particles with symmetric +/— unit charges. Thermal equilibrium is treated in the 
grand canonical ensemble characterized by the inverse temperature f3 and by the couple 
of equivalent (there is no external electrostatic potential), possibly position-dependent, 
particle fugacities z + (r) = z_(r) = z(r). In the disc geometry presented in figure 1, 
there are two domains: the disc of radius R, Ar = {r, |r| < R}, and its complement 
Ajj = {r, |r| > R}. The electrolyte in is modelled by a two-component plasma with 
the particle fugacity denoted by z. The wall A R is modelled by another two-component 
plasma with the particle fugacity z^, the choice z > corresponds to a conducting wall 
(e — > oo) while z = corresponds to vacuum (plain hard wall with e = 1). In 2D, the 
Coulomb potential v at a spatial position r, induced by a unit charge at the origin 0, is 
the solution of the 2D Poisson equation 

Av(r) = -2n5(r). (1.3) 

Explicitly, 

v(r) = -ln(|r|/o) (1.4) 

where a is a free length scale, which only determines the zero of the potential and should 
not enter statistical mean values. The interaction energy of charged particles {i,qi}, 
immersed in a homogeneous medium of dielectric constant = 1, is J2i<j QiQj v (\ r i ~ T j\)- 
The paper is organized as follows. 

In general section 2, the Coulomb gas confined to some 2D regions is shown to be 
equivalent, in the ideal-conductor limit, to the massless Gaussian model defined in the 
complementary empty regions. Since in the derivation |T3] of the large-size expansion of 
the free energy (|l.lj) for the critical Gaussian model with c = 1 only the curvature of the 
boundary is used, this explains the difference in sign of the logarithmic term between 
ijO and (fOl. 

In section 3, the underlying 2D fully-finite Coulomb system is solved in the Debye- 
Hiickel limit. All possible realizations of the model are considered, and the cancellation 
of the universal Casimir terms is documented when both the disc domain and its 
complement are occupied by a Coulomb gas. Fundamental differences between fully- 
finite and semi-infinite geometries are pointed out. 
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The exact solution of the model system at the Thirring free-fermion point (3 = 2 is 
presented in section 4. It is shown that basic features of Casimir phenomena, predicted 
by the mean-field theory, persist also at this specific temperature. 

A brief recapitulation is given in section 5. 



2. The Gaussian model approach 

The universal term (x/6) lni? in the grand potential (jl.2J) of a finite 2D Coulomb system 
closely resembles the logarithmic term in (jl.lj) valid for critical systems. At first sight, 
it is surprising that a Coulomb system, with short-range particle correlations, exhibits a 
critical-like behaviour. It has been argued [TBI E! that this is related to the existence of 
long-ranged critical-like correlations of the electric potential; reasons for the difference 
in sign of the logarithmic term between (jl.lj) and (jl.2j) have been given. A slightly 
different, and we believe more convincing, argument will now be presented. 



2.1. General formalism 

In the present subsection, we consider a 2D manifold on which some regions C are 
classical ideal conductors and some regions E are empty. The ideal conductors are 
obtained as the infinite fugacity limit of conductors having some microscopic structure, 
which amounts to take the limit in which the microscopic length scale goes to zero. 
Thus, the results which will be obtained in this limit are expected to be valid when 
the length scales under consideration are much larger than the microscopic lengths. 
For simplicity, we restrict ourselves to an infinite-plane manifold, CUE = R 2 , and 
start with conductors which are symmetric Coulomb gases of the kind described in the 
Introduction, but some generalizations are straightforward. 

The interaction energy E^ of ±1 charged particles {qj, ^j}jL\ is expressible in terms 
of the microscopic charge density p7v(r) = J2j=i Qj$( r ~ r j) as follows 

E N ({ qj , rj }) = ~ J r2 d 2 r J r2 d 2 r'p N (r)v(\v - r'\)p N (r>) - itf w (0), (2.1) 

where v (0) is the (diverging) self-energy. The thermodynamic properties of the Coulomb 
system are determined by the grand partition function H defined as follows 

oo oo N + N- 

2= E E ^J^Q(N + ,N.), (2.2a) 

N+=0N-=0 JV + - JV -- 



where 



Q(N + , N_) = / f[ d 2 r, exp \-f3E N {{q v r,-})] {2.2b) 

is the configuration integral of -/V+ positive and iV_ negative charges, iV = iV + + iV_ 
and z + = Z- = z are the equivalent fugacities of charged particles. Let us insert 
the energy representation (j2.1j) into (j2.2fcjl . The self-energy simply renormalizes the 
particle fugacities. With regard to the fact that, according to (jl.3j) . — A/(27r) is the 
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inverse operator of the Coulomb potential v(r), the standard Hubbard-Stratonovich 
transformation (see e.g. ref. [2E]) can be used to express the integral bilinear term as 

I a r 

/X?0exp [f R2 d 2 r (poAo + \ 



exp 



2 Jb? ^ Jb? d2r 'pN(r)v(\r - r'\)p N (r') 



N 



(2.3) 



/P0exp(/ fi2 d 2 r£0A0) 

where 0(r) is a real scalar field and / V<p denotes the functional integration over this 
field. The consequent factorization of the contributions from (N + , i\L) particle states 
in (|2.2a|) then allows to express the grand partition function of the system in the form 



/P0exp(-^(r)]) 
JV<f>eM-S[0}) 



with 



S[z(r)} 



d 2 r 



CUE 



Air 



0(r)(-A)0(r) - 2z(r) cos(/30(r)) 



being the 2D Euclidean action of the sine-Gordon theory and 

J z for r G C, 



z r 



for r e E. 



(2.4a) 



(2.46) 



(2.5) 



The normalization of z is fixed by the short- distance expansion of the two-point 
correlation function 



) 



-P 



as r 







(2.6) 



under which the self-energy factor disappears from statistical relations; for more details 
see ref. [30J and the references cited therein. 

Since the scalar 0-field has to be regular at infinity, the term </>(— A)</> in (|2.45jl 
can be transformed via the integration by parts into |V0| 2 . The sine-Gordon action 
(12.46j) thus takes its minimum at a 0(r) constant in space. Due to a discrete symmetry 
— * + 27m/(3 with any integer n, the action has infinitely many ground states |0 n ) 
characterized by the associate expectation values of the field ((f)) n = 27in//3. These 
ground states become all degenerate when the size of the Coulomb domain \C\ is large 
|3T] . It is therefore sufficient to develop the sine-Gordon action (|2.46jl . on the classical 
level as well as on higher quantum-correction levels, around any one of these ground 
states say |0 ) with (0) o = 0. Now, in the ideal-conductor limit z — > oo, the 
minimum of — 2zcos(/3(j)) around = becomes infinitely sharp and 0(r) identically 
vanishes in the regions C. Up to an irrelevant multiplicative constant, the numerator of 
(|2.4a|) becomes 

" A 

47T 

i.e. the partition function of the massless Gaussian model in the empty region(s). The 
boundary condition for is that it vanishes at the interface (s) of E and C, since it 



Zq = / V(f) exp 



d 2 r0(r)(-A)0(r) 



;2.7) 



Screening of universal Casimir terms 



7 



has to be continuous and vanishes in C. The denominator of (|2.4gj) is just an (infinite) 
constant, independent of the geometry of the regions E and C. 

Although has some resemblance with the electric potential, it also has some 
drastically different properties. For instance, in a bulk Coulomb gas, the of the sine- 
Gordon representation has short-range two-point correlations while the electric potential 
has long-range ones [32] ■ Here, the electric potential has fluctuations in the conductor(s) 
(these fluctuations survive in the ideal-conductor limit) and by continuity the electric 
potential has fluctuations at the interface(s) of E and C, while is identically zero at 
the interface (s). 

2.2. Parallel conducting plates in v dimensions 

In the case of two parallel conducting plates, at a distance d of each other, separated by 
vacuum, the Gaussian model ([2.7)1 . generalized to v dimensions, has been shown |1U| ITT] 
to have in its free energy per unit area of one plate a universal <i-dependent term F such 
that the corresponding electrostatic Casimir force per unit area —dF/dd is attractive 
and given by 



where s v is the area of the unit sphere in v dimensions and where T and ( are the gamma 
function and the Riemann zeta function, respectively. In the Appendix, we give another 
derivation of (|2.8|) for the Gaussian model. In subsection 2.2. of I, we have obtained the 
same Casimir force ()2.8[) by using the fluctuations of the electric potential (which do 
not vanish at the surface of the plates) . This is a check about the validity of the present 
method using the Gaussian model with simple Dirichlet boundary conditions for the </> 
field. 

2.3. A hole in a 2D ideal conductor 

We now come back to two dimensions. We consider an ideal conductor filling the whole 
plane except for an empty hole of characteristic size R, with a smooth boundary. It has 
been shown ^3] that the Gaussian model has a free energy of the form ([Lip with c = 1 
and % — \. Thus, (3 times the grand potential of the system has the universal term 
— (1/6) \nR. If the hole is a disc in a plane, R may be chosen as its radius. 

2.4- A 2D Coulomb gas surrounded by an ideal conductor 

We now consider the opposite geometry of an ideal conductor filling some connected 
region C of the plane, with a smooth boundary, surrounded by vacuum. Now the 
relevant Gaussian model fills the exterior of C. In the derivation of ref. [Tl], the 
curvature of the boundary is used, and it is easy to see that the present geometry differs 
from the one of Section 2.3 by a change of sign of the curvature, leading to the universal 
term +(1/6) mi? in (Oil . 




(2-K) u - 1 (2dY 



(2.8) 
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3. Debye-Hiickel theory 

3.1. General result 

The high-temperature limit of the model system in figure 1 is described by the Debye- 
Hiickel theory For the general formalism, see e.g. section 3.1. of paper I. The formalism 
applied to the 2D symmetric Coulomb gas can be briefly summarized as follows. 

In the bulk (homogeneous) regime, the total particle number density n as a function 
of the fugacity z is given by 

n l-(/?/4) 



za 



0/2 




(3.1) 



where C is the Euler's constant. The inverse Debye length is defined by k 2 = 2irj3n. 

In the inhomogeneous regime, the whole system domain A can be separated into 
disjunct physically non-equivalent subdomains, A = U Q A^ Q ^. Within the grand canonical 
formalism, each subdomain is characterized by a constant fugacity, z(v) = z^ for 
r G A( Q ); the choice = corresponds to a vacuum subdomain with no particles 
allowed to occupy the space. The corresponding "bulk" particle density n^ a ' is related 
to the particle fugacity z^ a ' via (|3.1|) . and the corresponding inverse Debye length is 
n a = (27r/3?7,( a ^) 1 / 2 . One introduces the screened Coulomb potential G which obeys 
within each domain A^ the differential equation 

A: - 4] G(n, r 2 ) = -27r<f(ri - r 2 ), n G A w . (3.2) 

Here, the spatial position of the source point r 2 is arbitrary. These equations are 
supplemented by the usual electrostatic conditions at each subdomain boundary dA^ : 
G and its normal derivative with respect to the boundary surface, d±G, are continuous 
at dA^ a \ The leading /^-correction to the constant particle density in the subdomain 
A( Q ) is then determined by linearizing the exponential in the expression 

n {a \r) =2z (a) exp{^ lim [-G(r, r') + v(|r - r'|)]l (3.3) 

2 r'^r J 

where the 2D Coulomb potential v is defined in equation (jl.4j) . 

We now apply the above formalism to the geometry of interest presented in figure 
1. The inverse Debye length will be denoted by k for the disc domain Ar = {r, r < R} 
and by Kq for the complement wall domain A^ = {r, r > R}. We first consider the case 
when both k and k,q are nonzero. Let the source point r 2 be first in the disc domain, 
i.e. r 2 < R. Equations ()3.2j) then take the form 

"Aj - k 2 ] G(r u r 2 ) = -2n5( ri - r 2 ), (r a < R), (3.4a) 

Ax - 4] G(n, r 2 ) = 0, (n > R). (3.46) 

In terms of polar coordinates, the solution of these equations can be written as an 
expansion of the form 

oo 

G(ri,r 2 )= J2 [^( Kr <)^( Kr >) +a l Ii(Kr 1 )Ii(Kr 2 ))exp[il(cp 1 - y? 2 )], (r x < R), (3.5a) 
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oo 

G(ri,r 2 )= bi K i( K on)Ii(Kr 2 ) exp[i/(y3i - ip 2 )], (n>R). (3.5b) 

Here, r < = min{r 1 ,r 2 } and r > = max{r 1 ,r 2 }. l\ and K\ are modified Bessel functions 
possessing the symmetry I; = if; = K_i. They satisfy the same differential equation 

f" + \f ~ (l + ^) / = 0, / = J, (re) or (3.6) 



but exhibit different asymptotic behaviours 

y/2lTX x-*oo \2x 



e x / 7T \ x / 2 

/,(*) ~ -z=, „ (t-) e"* (3.7) 



''<■'■> T7TT - St(T) • (3. 



and 

+ o|/|! V27 ' 1V 7 ^o2|/| 

except for the special Z = case of Kq(x) ~ — ln(a;/2) — C in the limit x — > where C is 
Euler's constant. We see that, in representations ()3.5ajl and (J3.5fej) of G, these asymptotic 
behaviours ensure the regularity of G at the origin and at infinity. The coefficients a; 
and bi are determined by the above defined boundary conditions for G(ri, r 2 ) at r\ = R: 

K^kR) + aJ^R) = biK^KoR), (3.9a) 
k [K[(kR) + chI&kR)] = KobiK'^KoR). (3.9b) 
Using the recursion formulas for the modified Bessel functions 

xl[(x) = xli±i(x) ± lli(x), (3.10a) 
xK[(x) = -xK l±l (x) ± IKi(x), (3.105) 

and the Wronskian relation 



k(x)K l+1 (x) + J, +1 (x)tf,(a:) = (3.11) 

x 

equations ()3.9a|) and (|3.96j) give 

a ' = — FT — Fa" + — — — — — , (3.12a) 



k(KR) RWi k(nR) 
1 



(3.126) 



RWi 7 

where the auxiliary quantity Wi is given by 
Wi = kIKkR^kvR) - Ii (kR) KqK[(kqR) 

= kIi +1 {kR)Ki(koR) + Ii(kR)k K 1+1 (k R) (3.13) 

= kI,_i (kR)Ki(kqR) + Ii(kR)koK^i(koR). 

Note the symmetries a\ = a_/, 6/ = 6_; and Wi = W-i. With regard to the differential 
equation (|3.fi|) obeyed by the modified Bessel functions, it is easy to show that Wi fulfills 
the equality 

^ln (RWi) = ^ 2 - kDUkR^KoR). (3.14) 
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The same procedure can be applied when the source point r 2 lies outside of the disc, 
i.e. r 2 > R. For r\ > R, one obtains the solution of the form 



G(ri,r 2 )= [h^or^Ki^Kory) + ciKi(K r 1 )Ki(K r 2 )]exp[il(ip 1 - ip 2 )]. (3.15) 

1= — oo 

The coefficients q are determined by the boundary conditions for G at r\ = R as follows 
c, = + 4,4^4. (3.16) 



A',(K„fi) BW,Jf,(KoH)' 

They fulfill the symmetry q = c_;. 

To obtain the density profile, note that in relations ()3.5aj) and ()3.15|) the first 
terms in the sums over I correspond to the expansion of the modified Bessel functions 
K (n\ri — r 2 |) and K (K \ri — r 2 |), respectively. In formula ()3.3|) for the particle density, 
these bulk contributions, minus the pure Coulomb potential v, imply the density- fugacity 
relationship (JHHJ). The second terms in the sums over I in equations (j3.5oj) and ()3.15j) 
are "reflected" contributions due to the boundary at r = R. After the linearization of 
(13.3)1 in (3, they lead to 



n(r) =n - —J2^ a i l i 



1=0 



n(r) = no - ^ ^fuciK? («or) 



(r < 
(r > #). 



(3.17a) 
(3.175) 



1=0 



Here, /x; is the Neumann factor: //q = 1 and [ii = 2 for / > 1. 

Having at one's disposal the density profile, the grand canonical partition function 
Sr(k, Kq) can be deduced in the following way. For the present geometry, S# is defined 
by 

N 



-R 



E 



i 



n 



NJN ' 



R 



d r j z + 



^ d 2 rj z ^J exp 



i<j 



(3.18) 



where iV + (iV_) is the number of positively (negatively) charged particles and iV = 
iV + + iV_. The averaged particle density at position r is given by n(r) = <K r — r i))- 
Combining this with equation (|3.18|) . one finds that 



n(R-) - n(i? H 



(3.19) 



where R~ {R + ) means the left (right) limit to R. With respect to the density profile 
[equations (|3.17qJ) and ()3.176|) with coefficients ai and q given by (j3.12aj) and (|3.16l) . 
respectively], this relation takes the explicit form 

d _ 1 00 1 

— lnS fl = 2nR(n - n ) - -{n 2 - kI)^2/j,i—Ii(kR)Ki(k R) 



Wi 



+fX> [k 2 Ii{kR)Ki{kR) - KqIi(koR)Ki(k R) 



(3.20) 



1=0 



The integration with respect to R of the second term on the rhs of ()3.20j) can be done by 
using ()3.14|) and the integration of the last term follows from the indefinite integration 
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formula for Bessel functions [3"i| 

2 

/ dxxIi(x)Ki(x) = ^- [2Ii(x)Ki(x) + I l+1 {x)Ki^{x) + Ii^{x)K l+1 {x)} 

= 1 -[(x 2 + l 2 )I l (x)K l (x)-x 2 r i (x)K' l (x) 
Thence, from ()3.20|) one gets 

In E R (n, n ) = const + -kR 2 (u 

1 oo 

+7l> {[M) 2 + /2 ] Ii{KR)Ki{KR) - (kR) 2 I[(kR)K[(kR)} 



(3.21) 



1 oo 

-5>ln(jnV,) 



(3.22) 



i=0 

oo 



1 oo 

7 X> { [(^) 2 + IiiKoR^KoR) - (k R) 2 I[(k R)K[(k R)} . 
4 z=o 



The integration constant is fixed by considering the R — > limit. With the aid of the 
asymptotic formulas ()3.8j) one gets, for instance, \im R ^ (RWi) = (k/ko)^. After simple 
algebra, the final result reads 



RWA^ 



In E R (n, no) = In E R=0 + nR 2 (n - n ) - - ^ /^ ln 

1 1=0 

-I oo 

{[(kR) 2 + I 2 ] Ii{kR)Ki{kR) - (kRYiKkR^KkR) - l) (3.23) 

4 ;=o 

1 oo 

--Ew{ [( K o#) 2 + ^] Iii^K^R) - {KoRfl'^^K'^R) - /} . 
4 z=o 

Here, E r= q is the grand canonical partition function of the system with zero disc radius, 
given by the obvious relation limA-K>o hiS# =0 /|A| = /3p(n ) = n [l — ((3/4)]. Formula 
()3.23|) was derived under the assumption that both particle densities n and n are 
nonzero. 

The limit n — > corresponds to no particles present inside the disc, i.e. the vacuum 
disc hole surrounded by the fluctuating conductor wall. Taking the n — > limit in ()3.23|) . 

one gets 



lnE fl (n = 0,n ) 



oo 

In E R=0 - nR 2 n - - ^ /ijn 
1 1=0 



(k q R) 



(3.24) 



1 oo 

-7 E VI { [(^) 2 + / 2 ] IiiKoR^KoR) ~ (KoRflliKoR^iKoR) - /} . 

4 z=o 

The limit n — > corresponds to no particles present outside of the disc, i.e. the 
Coulomb system in the disc surrounded by vacuum (the plain hard wall). One can take 
the n — > limit in every term of equation (|3.23|) . except for the I = term in the first 
summation on the rhs of that equation. This I = term makes a problem due to the 
logarithmic divergence of Kq(k,qR) in Wq in the limit — > 0. The problem with the 
I = term, caused by the fact that the effective Coulomb potential G(r) is not screened 
at asymptotically large distances r in the present limit k q — > 0, was discussed e.g. in 
|33]). According to this reference, when calculating the particle density (j3.17aj) one finds 
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that the coefficient a depends on the free length scale a of the 2D Coulomb potential 
(11.4)1 . Since a statistical mean value should not depend on a, the limit a — > oo, which 
puts the zero of the Coulomb potential to infinity, has to be considered. Our formalism 
in the limit k — > is equivalent to that outlined in if one sets k = 1/a. Thus, 
returning to the particle density (|3.17fl| . the auxiliary quantity W\ ()3.13j) . considered in 
the limit k,q — > and for R > 0, behaves like 

kI 1 _ 1 (kR)K 1 (k R) for / > 0, 



W, 



and the coefficients ai, given by ()3.12aj) . take the form 
K^(kR) 



a, 



/ = 0, 1, 



Ii-i(kR) ' 

On the basis of (|3.19)l taken with n(R + ) = 0, it holds 



d 1 - < 

—In n R {n,n 



0) = 2-nRn - 



K 2 i?~ K^kR) 



Vi-j— 



(3.25) 



(3.26) 



(3.27) 



The integration of this equation finally implies 



lnS^(n, n = 0) = 7ri? 2 n In + const 



oo 



(i-iy.i^R) 



(3.28) 



-1 oo 

[(kR) 2 + I 2 ] Ii(kR)Ki(kR) - (kR) 2 I[(kR)K[(kR) - l] . 

4 1=0 



Here, since the grand partition function depends on the length scale a = the 
integration constant is in fact infinite in the considered limit a — > oo. It is clear from 
the derivation procedure that formula (|3.28)1 is valid only for R > 0, and it cannot serve 
as a basis for an expansion around the R = point. 



3.2. Large-R analysis 

Each of the above derived grand potentials Q = — (l//3)lnH is given in terms of infinite 
sums which cannot be summed up explicitly. What can be done is the evaluation of 
first few terms of the asymptotic expansion of the sums for large disc radius R — > oo. 
Denoting by a either of the dimensionless combinations kR or k,qR, one has to use 
the Debye expansion [36J of the modified Bessel functions Ii(a),Ki(a) and of their 
derivatives, since this expansion is valid for large I uniformly with respect to a/ 1. In 
particular, one has 



Ida) 



K,(a) 



2t(« 2 + / 2 )!/4 
1 



1 + 



3t - 5t 3 



a- 



+ / 2 )V4 



24/ 
3t - 5t 3 

24:1 



+ 



a- 



+ l 2 



O 



a" 



P 



(3.29a) 
(3.296) 
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and 



Ilia) 



1 ( a 2 + f 2)l/4 

/2tt 



a 



-9t + 7i 3 



K[{a) = - 



a 



1 - 



24/ 

-9t + 7t 3 
247 



+ 



1 



+ 



a 2 + l 2 
1 



a 2 + / 2 



(3.30a) 
(3.306) 



where 



Tl(l, a) = Va 2 + l 2 -l sinh" 1 (I /a) , t = l/y/a 2 + l 2 . 



The consequent sums over / can be performed by applying the (generalized) Euler- 
MacLaurin summation formula 



E/(0 



, /(/)dZ+i[/(n)+/(m)]+^[/ , (n)-/'(m)]- 



£4 
4! 



-/"»]+■■ -,(3.31) 



where 5^ are Bernoulli numbers: B 2 = 1/6, B4 = —1/30, etc. 

In the case of the Coulomb gas of particle density n , localized outside of the disc 
of radius R with vacuum in the disc hole, formula 1)3.24)1 evaluated in the R — > 00 limit 
implies 

(3tt R (n = 0, n Q ) = -Pp{n ){\A\ - nR 2 ) + ^ 7 (n ) 2tcR - hn^R) + 0(1). (3.32) 

o 

The first term on the rhs of ()3.32)) is the bulk contribution with the pressure p given by 



Pp(n Q ) 



1 - 



4 



no, 



(3.33) 



the second term is the surface contribution with the surface tension 7 given, in the 
Debye-Hiickel limit, by j2Zj 



dl 

4:71 



hi 



(I 



8tt 



(4 



7T 



(3.34) 



4l^ 2 + I 2 

Finally, the third logarithmic term has the universal coefficient —1/6. This Casimir 
term, caused by electrostatic fluctuations inside the wall, tends to dilate the empty disc 
domain. This is the fundamental difference in comparison with semi-infinite geometries 
where fluctuating walls attract one another. 

In the case of the Coulomb gas of particle density n, localized inside the disc of 
radius R and surrounded by vacuum, formula (j3.28j) evaluated in the R — > 00 limit gives 



pn R (n, n = 0) = -Mn) txR 2 + ^7(71) lnR + -InUR) + 0(1) 

6 



(3.35) 

As before, the first and second terms on the rhs of 1)3.35)1 correspond to the volume and 
surface parts of j3Q, respectively. The third logarithmic term has the coefficient 1/6, 
in agreement with the general relation (jl.2)l taken at the disc value of x = 1- This 
coefficient is opposite to the one in 1)3.32)1 which confirms the relation of the universal 
term to the curvature of the constraining surface. It is important to note that in 
semi-infinite geometries (see e.g. paper I) the grand potential of a Coulomb system 
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constrained by vacuum plain hard walls does not exhibit the universal Casimir term. 
This is another fundamental difference between fully-finite and semi-infinite geometries. 

From a technical point of view, it is useful to sum the two expressions for the grand 
potential ()3.24j) and (|3.28|) . taken at the same particle density denoted say by ii: 

l3VL R {n = n, n = 0) + /3fl R (n = 0, no = n) = /3Q R= o + const 

OO r 1 

[RRfKi (RR)h (RR)] 



4" 



21 



[R,R) 2 Ii-i(kR)Ki + i(kR) 



(3.36) 



Using first the recursion formulas (3.10) and subsequently the asymptotic expansions 
(3.29) and (3.30) for the modified Bessel functions, the application of the Euler- 
MacLaurin summation formula ()3.31|) implies after simple algebra that 

f3VL R (n = n, n = 0) + (3Q R {n = 0, n = n) = f3Q R=0 + 2f3^(n)2-KR + 0(1). (3.37) 

This relation proves the consistency of asymptotic expansions ()3.32|) and (|3.35J) . and 
enables us to write down one knowing the explicit form of the other. Simplifying 
technicalities of this kind will be used in what follows. 

We are now ready for studying the large-i? asymptotic behaviour of the grand 
potential when both particle densities n and no are nonzero, see formula ()3.23|) . One 
has for the specific combination of grand potentials 

1 



/3Sl R (n, n ) - \J3n R (n, n = 0) + (3n R (n = 0, n )] 

IkRI^kR^oRK^oR)} 



const H — In 
2 




W 

kk RIi(kR)K 1 (k R) 

aw, j ■ (3 38) 

As above, using the recursion formulas and the asymptotic expansions for the modified 
Bessel functions, the Euler-MacLaurin summation formula leads to 

f3tl R (n, no) - lPQ R (n, n = 0) + (3Vt R {n = 0, n )] 

21 (^MTP + M + I 2 



R 



dl In 



0(1). 



(3.39) 



(/ + Vk^TP) (l + yjK 2 + 1 2 
With regard to the asymptotic expansions ()3.32|) and ()3.35|) . one finally arrives at 
(3n R (n,n ) = -{3p(n)irR 2 - {3p(n )(\A\ - nR 2 ) + (3^(n,n )2nR + 0(1), (3.40) 
where 



00 dl, 
— In 
An 



A^/k^+Pjnl + 1 



(3.41) 



is the (dimensionless) contact surface tension of the two 2D plasmas in the Debye-Htickel 
limit, see equation (3.26) of paper I. The universal logarithmic Casimir term does not 
appear in ()3.40j) . it is "screened". This phenomenon is intuitively expected: since there 
are Coulomb systems on both sides of the disc boundary, the curvature contributions 
with opposite signs cancel with one another. The same cancellation of long-ranged 
Casimir forces takes place in semi-infinite geometries, see paper I. 
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4. The free-fermion point 

4-1. General result 

The 2D Coulomb gas of symmetric unit charges is exactly solvable at the collapse point 
(3 = 2; for general formalism see e.g. section 4.1. of paper I. In the grand canonical 
formalism, at /3 = 2, the bulk system is characterized by the rescaled particle fugacity 
m = 2naz [a is a free length scale introduced in p.2jl ] which has the dimension of an 
inverse length. The many-particle densities can be expressed in terms of specific Green 
functions G qq >(r, r') (q, q' = ±); because of the symmetry between positive and negative 
particles one only needs G ++ and These Green functions are determined by the 

equations 

(Ai - m 2 )G ++ {r 1 , r 2 ) = -mS^ - r 2 ) (4.1) 

and 

G_ + (n, r 2 ) = -- ( A + i A ] G ++ (n, r 2 ), (4.2) 



supplemented by the vanishing boundary conditions when |r x — r 2 | — > 00. In infinite 

space, the solution of (|4.1|) reads 

777 

G ++ (T 1 ,T 2 ) = —K {m\T 1 -T 2 \). (4.3) 

The one-particle densities n + = n_ = n/2 (n is the total particle density), given by 

n«W=mG M (r,r), (4.4) 

are infinite since K (mr) diverges logarithmically as r — > 0. Regularization of the 
Coulomb interaction by a short-distance cutoff L implies for the particle density 



m T . , _ * ml", / 2 
n = — K {mL) ~ — 

7T mL->0 7T 



In ( — ) - C 

mL 



(4.5) 



In the inhomogeneous case when the system domain A = U a A^, each subdomain 
A(°) is characterized by a constant rescaled fugacity, m(r) = m a for r G A^, and 
the corresponding bulk density n a defined as a function of m a by (|4.5|) . Within each 
domain, the Green function G++ obeys the differential equation 

(A 1 - m*)G ++ (ri, r 2 ) = -m Q 5(n - r 2 ), n e A^, (4.6) 

where the spatial position of the source point r 2 is arbitrary. G_ + is determined by 
relation ()4.2|) with m substituted by the sub domain-dependent m(ri). The boundary 
conditions are that both C++ and G-+ must be continuous at each subdomain boundary 
dA( a \ The one-particle densities are again given by (|4.4|) with m — > m(r). 

For the disc geometry of figure 1, the rescaled particle fugacity is equal to m in the 
disc domain r < R and to m in the complementary wall r > R. Let the source point 
r 2 be first in the disc domain, i.e. r 2 < R. Equations (j4.1j) then take the form 

(Ai - m 2 )G++(r 1 , r 2 ) = -m5(n - r 2 ), (n < (4.7a) 
(Ai - m^G ++ ( ri) r 2 ) = 0, (n > i?). (4.76) 
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In terms of polar coordinates, the solution of these equations is written as an expansion 
of the form 

m 00 

G++(ri,r 2 ) = — [H mr <) K l( mr >) 

Z7T , 

(=— 00 

+ail l {mr 1 )li(mr 2 )} exp[iZ(<pi - (p 2 )), (n < R), (4.8a) 

m 00 

G++(ri,r 2 ) = — biK l {m r 1 )Ii(mr 2 )exp[il{ip 1 - <p 2 )], {n>R). (4.86) 

271 , 

The polar version of relation (|4.2|) for G_ + takes the form 

e ifi ( Q l d \ 

GL + (r 1; r 2 ) = — — + - — G ++ (n, r 2 ). (4.9) 

myci) \ar\ r x o<fi J 

The coefficients a\ and bi are determined by the continuity conditions for G ++ (ri,r 2 ) 
and Cr_ + (ri, r 2 ) at the disc boundary r\ = R: 

Ki(mR) + ailiimR) = biKi(m R), (4.10a) 
K l+1 (mR) - aj l+1 (mR) = b t K l+1 (m R) . (4. 106) 

The solution reads 

KijmR) 1 K t (m R) 
ai Ii{mR) mRVi k(mR) ' ' 

b > = sk- (4114) 

where 

Vi = I l (mR)K l+1 (m R) + I l+1 (mR)K l (m R). (4.12) 

Note the symmetry V\ = VLj_i. It can be readily shown that the auxiliary quantity VJ 
fulfills the equation 

Q <ry) — Tft n 

—In (RVt) = [I^mR^moR) + I l+1 (mR)K l+1 (m R)} . (4.13) 

OK V\ 

One proceeds analogously when the source point r 2 lies in the wall, i.e. r 2 > R. For the 

case r\ > R, one gets 



G+ + (ri,r 2 ) = — - [h{.m r < )K l (m Q r > ) 

2TX , — 

(=— 00 

+ciK l (m r 1 )Ki(m r 2 )) exp[i/(y?i - cp 2 )), (n, r 2 > R). (4.14) 

The coefficients q are given by the boundary conditions at r\ — R as follows 

h(m R) 1 IjjmR) 

Cl Ki(moR) moRVi KiimoR)' 1 ' ] 

The density at r — > determined by (|4.4jl . (|4.8aj) and (|4.11qjl . reads 
m 00 1 

n(iT) = — E jrUmRlKtimoR). (4.16) 

(=— 00 ' 
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The sum in (J4.16|) is divergent, and so it has to be formally regularized by taking an 
upper cutoff on I. Similarly, the density at r — * R + , determined by (|4.4j) . (|4.14j) and 
(I4.15J) . is expressible as follows 



n (R+) = ^n(R-). 
m 



(4.17) 



The grand canonical partition function H#(m, mo) is again determined by the differential 
equation (|3.19jl . in particular 



d 00 1 

lnS^ = 2(m - m ) £ - [I l {mR)K l {m R) + I l+1 {mR)K l+1 {m R)) . 



dR 



Due to the equality (j4.13j) . this equation can be integrated to the form 

00 

In E R (m, mo) = const + 2 In(RVi). 



(4.18) 



(4.19) 



1=0 



The integration constant is fixed by the R — > limit. Since lim^o RVi — m o 1 (m/mo) z 
for I > 0, one gets 

lnS R (m,m ) = lnS R=0 + 2j^ln<| [ — J m R 



1=0 



m J 



x [I l (mR)K l+1 (m R) + I l+1 (mR)K l (m R)] \. 



(4.20) 



Here, S# = o is the grand canonical partition function of the system with zero disc radius, 
i.e. limA-KxJnHi^o/IAI = f3p(n ). 

Although (|4.20|) was derived under the assumption that both rescaled fugacities m 
and m are nonzero, there is no problem to consider the zero limit of either of fugacities 
directly in (|4.20|) . One obtains 



In E R (m, m = 0) = 2 ^ In 

1=0 



and 



In E R (m = 0, m ) = lnS^o + 2 ^ In 

1=0 



2 fm R\ l+1 



. (4.22) 



4-2. Large-R analysis 



The result (|4.21|) for the 2D two-component plasma at (3 = 2, in the disc surrounded 
by vacuum, has already been derived in jlbj . There, the large- R asymptotic form of the 
grand potential was derived in the form 

0) = -(3p(m)7rR 2 + (3j{m)2nR + -\n(mR) + 0(1), (4.23) 

6 

where p(m) is the regularized pressure and the surface tension 7(772) is given by 

n / \ ( 2Vm 2 + / 2 \ (\ \\ fAnA . 

= —In - ===== = m - - — . (4.24) 
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One can apply the standard procedure to formula (J4.22j) corresponding to the 
plasma outside of the empty disc, to obtain the large- R asymptotic behaviour 

(3Q R (m = 0, mo) = -/3p(m )(\A\ - irR 2 ) + f3-f(m )2nR - hn{m R) + 0(1). (4.25) 

o 

One sees that, similarly as in the Debye-Hiickel limit, in comparison with ()4.23|) the 
universal logarithmic term has the opposite sign. 

Technically, it is simpler to derive (|4.25jl by summing up the two basic expressions 
(I4.21|) and (|4.22|) taken at the same rescaled fugacity, say to: 

(3VL R {m = m, m = 0) + f3Q R (m = 0, m = to) 

oo 

= pn R=0 - 2^1n [mRIi(fhR)K l+1 (ffiR)] . (4.26) 

1=0 

The standard asymptotic procedure for the modified Bessel functions leads to the 
following asymptotic behaviour 



(3VL R {m = m, m = 0) + l3Vt R (m = 0, m = fh) = (3Q, R=0 + 2/?7(m)2vr J R + 0(1), 

which proves the consistency of relations (|4.23J) and ()4.25|) . 

Finally, considering the special combination of grand potentials 

(3Vt R (m, m ) - [flQ R (m, m = 0) + (3VL R (m = 0, to )] 

I l+1 (mR)K l (m RY 



(4.27) 



2Eln 

1=0 



1 + 



I l (mR)K l+1 (m R)_ 
and using the standard asymptotic procedure, one gets 
(3Q R {m, too) - \ftfl R (m, m = 0) + (3Vt R (m = 0, to )] 



(4.28) 



R 



dl In 



(l + Vm 2 + I 2 ) (l + ^/m 2 + l 



+ 0(1). 



2 ^mm + y/m 2 + l 2 ^fm 2 + I 2 + I 2 

With regard to the asymptotic relations (|4.23|) and (|4.25|) . one concludes that 
f3il R (m, m ) = -pp(m)nR 2 - pp(m ){\A\ - ttR 2 ) + ^(m, m )2nR + 0(1), 
where 



2Vm 2 + I 2 J ml + I 



mm 



+ Vm 2 + l 2 Jm 2 + l 2 + l 



(4.29) 



(4.30) 



(4.31) 



defines the surface tension of the two 2D plasmas at (3 — 2, see equation (4.12) of paper 
I. As was expected, the universal logarithmic term disappears once again. 

We have verified that the basic features of the mean-field behaviour, predicted by 
the Debye-Hiickel analysis in the previous section, persist also at the specific inverse 
temperature [3 = 2. 
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5. Conclusion 

The name "ideal-conductor" usually means that the electric potential has some constant 
value (for instance zero) inside the conductor, without any fluctuations. More precisely, 
we call such a model of a conductor "inert ideal-conductor". The present paper deals 
with microscopic models of conductors, in which there are fluctuations of the electric 
potential. They are "living" conductors. In the limit when the microscopic lengths, 
such as the Debye length, go to zero (in practice, are small compared to the macroscopic 
length, here such as the radius R), the fluctuations of the electric potential survive. We 
call a conductor in that high-density limit "living ideal-conductor" . 

In a previous publication J7|, Coulomb systems with inert ideal-conductor 
boundary conditions were studiedf. In this model, for a Coulomb system in a slab 
of width d, a repulsive Casimir force was found, which is just the opposite of (|2.8jl which 
holds for living ideal-conductor plates separated by vacuum. For a Coulomb system in 
a slab with living ideal-conductor walls, the two contributions cancel each other, leaving 
only a short-range attraction, as shown in I. 

Similarly, in ref. jl^j, a 2D Coulomb system in a disc of radius R with inert ideal- 
conductor walls was found to have a logarithmic universal contribution (1/6) InR to (3 
times its grand potential. This is just the opposite of the —(1/6) In R found here for an 
empty circular hole surrounded by a Coulomb gas. Again, for a Coulomb system in a 
disc with living conductor walls, the two contributions cancel each other. 
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Appendix 

In ref. [T7|, ()2.8|) has been derived by using the quantum Hamiltonian in v— 1 dimensions. 
Here, we give an alternative, more direct, derivation. 

The free energy per unit area F corresponding to the Gaussian partition function 
(12. 7|) . generalized to v dimensions, is given, up to an irrelevant additive constant, by 

^ = -7^T ln ^ = ^TElnA, (A.l) 

where L is the linear size of a plate and A the eigenvalues of minus the Laplacian. The 
planes, on which <fi obeys Dirichlet boundary conditions, are perpendicular to the x axis 
at x = and x = d. Each point is defined by its Cartesian coordinates (x, r ) where r 1 - 
is a {y — l)-dimensional vector normal to the x axis. In this geometry, in the limit of L 
infinite, the eigenfunctions of —A are exp(— il • r 1 -) s\.n(mtx/ d) where 1 is a wavevector 

\ In JJji the Coulomb potential in v dimensions {y > 2) was defined as r 2 ~ v while in I it was defined 
as r 2 ~ v ' j{y — 2). These different definitions do not change the Casimir terms under consideration. 
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normal to the x axis and n an integer larger than 0. The corresponding eigenvalues are 
A = I 2 + (nir/d) 2 . Therefore, (jA.ljl can be written as 



1 r d"' 1 ] °° 
(3F = - ^— l - Y In 

2 J (2tt) 1 '- 1 



/ 2 + 



n=l 




Thus 



P dd 



d y -H 1 ~ 



(mr/d) 2 



(A.2) 



(A.3) 



(2 7 r) I '- 1 d^ 1 I 2 + (mr/d) 2 ' 

The sum on n in ()A.3|) diverges. The divergent part can be separated by noting 
that the summand minus 1 is —(ld/7r) 2 /[(ld/Tr) 2 + n 2 ], which can be explicitly summed 
[HS]- The result for the total sum, including the divergent part, is 



1 ~ (nn/d) 2 
d^P + inir/d) 2 



--lcoth(ld) + — 
2 K J 2d 



OO 1 

n= 



1i d' 



(A.4) 



The last two terms of (jA.4|) can be regrouped into the divergent sum S = 
(1/2) H^L_ 0O (l/(i). Since F has a bulk term of the form Ad, with A infinite in the present 
ideal conductor limit, for obtaining the finite-d Casimir force /, we must subtract from 
(|A.3j) its value for d infinite. For d infinite, the term nn/d in the eigenvalue A becomes 
the continuous wave number k and the divergent sum S becomes dk / (2tt) . Since, 
for d finite, the summand in S is a constant, S can also be written as the same integral 
on k, and after the subtraction the divergent term disappears. Finally, 



pi 



dd dd 
d v -H 



d=oo. 



(2vr) 1 



T /[coth(/rf) - 1] 



(27r)^- 1 (2d)^ 



(A.5) 



in agreement with (|2.8|) . 
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